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ABSTRACT 



Aims. We present high spatial resolution observations of the mid-infrared core of the dusty symbiotic system HM Sge**. 
Methods. The MIDI interferometer was used with the VLT Unit Telescopes and Auxiliary Telescopes providing baselines 
oriented from PA=42° to 105°. The MIDI visibilities are compared with the ones predicted in the frame of various 
spherical dust shells published in the literature involving single or double dusty shells intended to account for the 
influence of the hot White Dwarf. 

Results. The mid-IR environment is unresolved by a 8m telescope (resolution~300 mas) and the MIDI spectrum exhibits 
a level similar to the ISO spectra recorded 10 yr ago. The estimated Gaussian Half Width at Half Maximum of the 
shell of 7.8±1.3mas (12AU, assuming a distance of 1.5kpc) in the 8-9/im range, and 11.9±1.3mas (18AU) in the 11- 
12/im range, are much smaller than the angular separation between the Mira and the White Dwarf of 40 mas (60AU). 
The discrepancies between the HWHM at different angle orientations suggest an increasing level of asymmetry from 
13 to 8/j.m. The observations are surprisingly well fitted by the densest (optically thick in the N band) and smallest 
spherical model published in the literature based on the ISO data, although such a model does not account for the 
variations of near-IR photometry due to the Mira pulsation cycle suggesting a much smaller optical thickness. These 
observations also discard the two shells models, developed in an attempt to take into account the effect of the White 
Dwarf illumination onto the dusty wind of the Mira. These models are too extended, and lead to a level of asymmetry 
of the dusty environment tightly constrained by the MIDI visibilities. These observations show that a high rate of dust 
formation is occurring in the vicinity of the Mira which seems to be not highly perturbed by the hot companion. 

Key words. Techniques: interferometric; Techniques: high angular resolution; Stars: AGB and post-AGB; Stars: binaries: 
symbiotic; Stars: circumstellar matter; Stars: mass-loss 



1. Introduction 

HM Sge is a D-type (dust forming) sy mbiotic system tha t 
erupted as a symbiotic nova in 1975 (|Dokuchaeval H976T) . 
evolving from a 17th to a 11th magnitude star with a rich 
emission line spectrum. The cool component is a Mira star 
and the hot component is a White Dwarf (WD) which had 
escaped detection enshrouded in the dense Mira envelop 
before this dramatic event. 

HM Sge was observed in IR from the ground and 
also with the IRAS and ISO satellites. The observa- 
tions show a spectrum dominated by silicate dust whose 
flux evolved deeply in about 30yrs. Many papers inves- 
tigated the time variable dusty environment and tried 
to evaluate the impact of the temper ature and luminos 
ity evol u tion of the WD on the dust dDokuchaevaL 1976; 



Wenzell. 119761; ICiatti et all Il977t iDavidson et all 



1978; 



Ciatti et all 119781; iKwok fc Purton Ul979tlBrvan fe Kwok . 



19911 : lYudin et all Il994fc IBogdanov fc Taranovar ]2001 



Send offprint requests to: S. Sacuto 

* Based on observations made with the Very Large Telescope 
Interferometer at Paranal Observatory under programs 075. D- 
0484 and 077.D-0216 

** Reduced visibilities and differential phases are 
available in electronic form at the CDS via anony- 
mous ftp to cd sarc.u-strasbg.fr (130.79.128.5) or via 



http:/ /cdsweb. u-strasbg.fr/cgi-bin/qcat?J/A-|-A/ 



ISchild et all l200lh . The strong infrared excess detected 
shortly after the outburst indicates the presence of a rather 
dense dust envelope. The dust content of the system was 
severely affected by the nova outburst and it was soon ob- 
serve d that a process of dissipation o f the dust was on the 
way dTaranova fc Yudinl. fl982l I1983D. de spite some obscu- 
ration events (|Munari fc Whitelockl . ll989l) . The IRAS spec- 
trum obtained in 1983 is about twice the ISO spectra ob- 
tained in 1996-1997 and this flux difference was attributed 
to a long -term variation effect and no t to a beam size 
variation. IBogdanov fc Taranoval ((2001) attributed these 
changes to an increase of the concentration and decrease 
of the temperature of the dust close to the Mira star, as- 
sociated with the temperature increase and l uminosity de- 
crease of the White Dwarf. In the other hand, ISchmid et al.l 
(2000) find that the red giant spectrum has strongly bright- 
ened and find that the dust column in the line of sight has 
markedly decreased. It is difficult to assess whether this is 
a long-term evolution or an other event related to the pre- 
sumably patchy circu mstellar environmen t. This last hy- 
pothesis is favored by ISchmid et al.l (|2000D and it must be 
noticed that longer wavelength emission is less sensitive to 
the dumpiness of the e nvironment which probably im plies 
that the conc l usion s of IBogdanov fc Taranoval (|2001h and 
ISchmid et all (|2000l) are not necessarily incompatible. 

Classically, most of the studies evaluated the infrared 
excess by fitting model spectral energy distributions (SEDs) 
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to the instantaneous observed ones. The common model 
is a single dust shell, whose composition is dominated 
by silica tes (but some s t udies involve a contribution from 
carbon (jBrvan fc Kwold . Il99lh ) . However, this simple ge- 
ometry could be considered as a crude approximation of 
a putatively complex geometry and dust composition of 
the mid-IR emitting regions. As an attempt to overcome 
this problem, more complex models involving two shells 
were used: a cool one surrounding the Mira star and 
a second, more extended that sees mostly the WD flux 
(|Bogdanov fe Taranoval l200llSchild et alll200lh . The best 
argument for a two shell model is that the dust-rich area 
emission is affected by the combined effects of a ioniz- 
ing UV irradiation from the hot primary on the expand- 
ing dust shell of the red giant, and a compression wave 
of the colliding winds in the transition zone. The conclu- 
sions of such approaches appear in contradiction to each 
other: whereas ISchild et alj (|200ll ) found a better match of 
the data with a two-compone nt model of dust distribution, 
iBogdanov fe Taranoval (|200lh are lead to the opposite con- 
clusion and claim that current models with only a single 
Mira-type source are preferable. It is important to notice 
that both studies use the same ISO spectra as constraints 
and that they are able to provide satisfactory fits of the 
SED with very different model parameters. It is obvious 
that the single constraint of a SED is not able to fully con- 
strain a putatively complex dust environment whose com- 
position, geometry and even the luminosity of the heating 
sources is poorly known. An additional difficulty is that the 
mid-IR observations are to date scarce and secured only 
with relatively extended beams. Moreover, HMSge is not 
in an equilibrium state (Mira pulsations, wind interactions, 
dissipation/formation of obscuring dust,. . . ) and this pic- 
ture might well have evolved over the past few years since 
the ISO observations in 1996-1997. There is also a need 
for a SED taken with a much smaller beam than the ISO 
one which would allow to isolate the signature of the most 
recently formed dust. 

We present in this article the first mid-infrared in- 
terferometric observations of this system taken with 
the MIDI/VLTI instrument recorded in 2005/2006 
providing a spectrum update from a 8m telescope 
beam (FWHM^300 mas) and a 1.5m telescope beam 
(FWHM~1.5"), and measurements of spatially correlated 
flux in interferometric mode that provide high spatial res- 
olution (~20mas) and multi-spectral (A/AA ~ 30) infor- 
mation on the mid-IR source. The flux of VLT single dish 
spectrum is directly comparable to the the ISO ones (see 
Fig. [T|). This is a clear indication that, first, most of the 
ISO emission observed 8yr before probably originated from 
a spatially compact region and second, that the dust con- 
tent of the system seems not to have changed significantly. 
As a consequence, we assume in the following study that the 
spatial dist ributions expected f rom the different models de- 
velope d by ISchild et al.l (|2001l ) and IBogdanov fe Taranoval 
(|2001h can be checked against these high spatial resolution 
observations and that the MIDI spectrum can be comple- 
mented with the ISO spectra to construct the current SED 
of the system at longer infrared wavelengths. The outline of 
the paper is the following: we present first the MIDI obser- 
vations in Sec (21 and we summari ze the published pa ram- 
eters of the models develope d by ISchild et al.1 (|2001h and 
IBogdanov fe Taranoval (|200lD in SecH Then in Sec|l the 
expected visibilities from the models are compared to the 



MIDI observations (the detail of the calculations is shown 
in appendix). We propose some improvements for the best 
model in Sec|5l We discuss the limitations of the approach 
based on spherically symmetric dust shells for the study of 
a spatially complex source in Sec [51 Finally we summarize 
and conclude in Sec[7l 



2. Observations 

The Very Large Telescope Interferometer (VLTI) of ESO's 
Paranal Observatory has been used with MIDI, th e MID - 
infrared Interferometric recombiner (jLeinert et all [2003). 
MIDI combines the light of two telescopes and provides 
spectrally resolved visibilities in the N band atmospheric 
window. 

The observations of HM Sge were conducted with the 
VLT Unit Telescopes (UTs) UT2, UT3 and UT4, providing 
projected baselines in the range of 32-59 meters and the 
Auxiliary Telescopes (ATs) E0 and GO, providing projected 
baselines in the range of 13.5-16 meters. 

The UTs observations were made during the nights of 
July 23-24th 2005, May 17th 2006 and June 11th 2006. The 
phase of the Mira during the UTs MIDI obs ervations was 
0/75 in 2005 and between 0.31-0.36 in 2006 (|Yudin et all 
1994). The data were recorded with different projected 
baselines, (46.5m, 44°), (32.1m, 42°), (37.8m, 47°), (46.2m, 
47°), (59.3, 105°) and (46.8, 101°). The ATs observations 
were made during the nights of May 21th 2006, May 27th 
2006 and June 16th 2006. The phase of the Mira during 
the ATs MID I obse rvations was between 0.32-0.37 in 2006 
(jYudin et all I1994D . The data were recorded with differ- 
ent projected baselines, (13.5m, 80°), (15.8m, 69°), (15.8m, 
75°) and (15.3m, 77°) but the night of May 21th 2006 was 
discarded due to the misalignment of the star (see below) . 
All the observations were made under good atmospheric 
conditions (seeing ~1 .01") with the worst seeing (^1.75") 
during the night of July 23th 200© 

The data reduction software^] MIA and EWS (|Jaffel 
12004ft were used to r educe the spectra and visibilities 
(jChesneau et all 120051 ). Chopped acquisition images are 
recorded (f=2Hz, 2000 frames, 4 ms per frame, 98 mas per 
pixel) for the fine acquisition of the target. The acquisi- 
tion filter was a N band filter. The prism of MIDI was used 
providing a spectral dispersion of about 30. 

Tab. Q] presents the journal of the interfero- 
metric observations. The calibrators, HD188512 
(G8IV, diam=1.98±0.02mas), HD187642 (A7V, 
diam=3.22±0.01mas) and HD206778 (K2Ib, 
diam=8.38±0.09mas) were observed right before or 
after each science target observations. 

The spectral template from ICohen et al.1 (|1999f ) of 
HD180711 (G9III) was scaled to the 12/xm flux of the in- 
terferometric calibrator HD188512 for absolute flux cali- 
bration. In interferometry, in parallel with measurements 
of the source, measurements of the nearby calibrator are 
taken. Because the difference in airmass did not exceed 0.1, 
no airmass correction was applied to the corrected spectra. 



The Adaptive Optics device MACAO dedicated to the VLTI 
provide diffraction limited images in N band for the UTs even 
for the worst condition encountered on this observation 

2 http://www.mpia-hd.mpR.de/MIDISOFT/ 
|http://www.strw. leidenuniv.nl/~nevec/MIDI/ 
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Table 1. Journal of observations: MIDI/UT2-UT3/UT3- 
UT4/E0-G0. The calibrators used to calibrate the visibili- 
ties are given below the science target. The phase of the 
Mira (<fiMira) during the observations is indicated. The 
length and position angle of the projected baseline are also 
indicated. 



Star 


UT date & Time 


<PMira 


Base[m] 


PA [deg] 


HM Sge 


2005-07-23 06:48:42 


0.75 


46.5 


44 


HD188512 


2005-07-23 07:12:10 








HM Sge 


2005-07-24 02:41:11 


0.75 


32.1 


42 


HD188512 


2005-07-24 02:58:55 








HM Sge 


2005-07-24 03:42:18 


0.75 


37.8 


47 


HD188512 


2005-07-24 04:03:14 








HM Sge 


2005-07-24 06:14:49 


0.75 


46.2 


47 


HD188512 


2005-07-24 06:37:58 
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HM Sge 


2006-05-17 09:04:02 


0.31 


59.3 


105 


HD187642 


2006-05-17 09:27:14 








HM Sge 


2006-06-11 09:05:09 


0.36 


46.8 


101 


HD187642 


2006-06-11 08:42:52 








HM Sge 


2006-05-27 08:17:39 


0.33 


15.8 


75 


HD206778 


2006-05-27 08:39:57 








HM Sge 


2006-06-16 06:28:01 


0.37 


15.3 


77 


HD206778 


2006-06-16 05:53:35 
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I i i i I i i , I , i i I i i i 

6 8 10 12 14 

Wavelength in micron 

Fig. 1. The solid line with error bars is the MIDI flux 
(Mira pulsation phase, (p=0.75) and the dashed and dashcd- 
dotted lines are the ISO/SWS spectra (upper one: 1996, 
^=0.65, bottom one: 1997, <p=1.08) 



The MIDI spectrum compared to the ISO spectra taken 
at two different phases of the Mira is shown in Fig. [TJ The 
level of flux received by MIDI is equivalent to the flux level 
of ISO spectrum, implying that most of the N flux is con- 
centrated within the 300 mas beam of a 8m telescope and 
also that the level of emission has probably not significantly 
changed since the ISO observations (1996-1997). We notice 
that the flux error bars of the photometry from UTs are be- 
low 8% level. The shape of the silicate feature in the MIDI 
spectrum is different from that of the ISO spectra: the flux 
level at 8/xm is close to the upper ISO spectrum whereas the 
flux longward 12/im is closer to the lower ISO spectrum. It 




o o r I o oo l l o 

7 8 910 1112 13 7 8 9101112 13 

wavelength (/im) wavelength (fJ.m) 



Fig. 2. MIDI visibilities (error bars) for the six projected 
baselines. The solid lines correspond to the uniform disk 
angular radius in mas and the dashed lines correspond to 
the FWHM in mas of the Gaussian distribution (both to 
be read from the scale on the right axis), computed from 
V at each wavelength. 

is difficult to state whether this is due to the signature of 
the dusty source at this particular pulsation phase or the 
consequence of a longer term evolution. 

We used all calibrators of each observing night to eval- 
uate the quality of the calibration. The data were reduced 
with both MIA and EWS. For the UTs data, the quality is 
uniformly good and we decided to apply a uniform multi- 
plicative error of 10% as a conservative estimate. The data 
reduction of the ATs data was much more difficult for sev- 
eral reasons: the alignment of the source was never optimal 
and strong biases in the extracted photometries were ob- 
served, affecting more the edges of the N band (7.5-9^m and 
12-13. 5^m). Moreover, the data reduction softwares are not 
currently optimized for these telescopes and developments 
are under way. As a consequence, we decided to show these 
data, but put only little weight on them compared to the 
more robust UTs data in the interpretation process. Some 
MIDI visibility curves were merged and averaged. We con- 
sider the two UTs observations of the July 23th 2005 at 
06:48 and July 24th 2005 at 06:14 as approximatively re- 
dundant because of their small baseline and position angle 
(PA) variation (<1% and <7%, respectively). This is also 
the case for the 2 ATs observations (May 27th 2006 and 
June 16th 2006) which have been averaged due to their 
small projected baselines and also their bad quality. The 
data from May 21th 2006 were disregarded due to alignment 
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problems. The other ATs data were reduced, calibrated and 
merged in a single curve in Fig. [21 We notice that the error 
bars at the edges of the N band, strongly affected by the 
photometric biaises are probably underevaluated. The six 
remaining calibrated visibility data sets of HM Sgc over the 
wavelength range from 7.5 to 13/im are shown in Fig. [2] 
The figure also shows the plot of the equivalent uniform- 
disk angular radius and FWHM of the Gaussian distribu- 
tion calculated from the visibilities at each spectral channel. 
The size of both models increases from 7.5 to 10/zm and is 
approximately constant between 10 and 13/im. Notice that 
the wavelength range has been slightly increased to show 
the MIDI visibilities between 7 and 8/im, at the edge of 
the atmospheric window. The quality of the data is good, 
and we can see a decrease of visibility that corresponds to 
the very strong emission of the forbidden line [NeVI] at 
7.65/mi (see Fig. Q]). This narrow line is observed in the 
best quality MIDI spectra but almost disappeared in the 
mean spectrum presented in Fig. [T]due to small mismatch 
of the source position on the slit. This decrease can be in- 
terpreted as a much larger extension for the line emitting 
region. This is not investigated further in the paper due 
to the limited spectral resolution and the potentially large 
errors at this wavelength. Finally, all visibility and differ- 
ential phase data as well as all the characteristics of the ob- 
servations are available from the CDS (Centre de Donnees 
astronomiques de Strasbourg) ; all data products are stored 
in the FITS-based, optical int erferometry data e xchange 
format (OI-FITS), described in lPauls et all (l2005f) . 

3. Spectrophotometric models 

The infrared observations (mostly from IRAS and ISO) 
were interpreted within the fra me of spherical mod els gener- 
ated by the DUSTY package (jlvezic et all 1 19991). We first 
reprod uce the SED s of models published b v ISchild et al.1 
(|2001h (SOI) and iBogdanov fc Taranoval (|2001h (B01). 
Then, still with the DUSTY code, we generate the cor- 
responding synthetic spectrally-dispersed visibility profiles 
throughout the N band (7.5-13/im) and compare them with 
the MIDI visibilities for each baselines. The aims of our 
study are threefold: first, finding the model that provides 
the closest match with the data; second improving, if pos- 
sible, the model parameters; and third, based on the best 
model, investigating a putative departure of spherical sym- 
metry of the source. 

3.1. Description of the models 

The models generated by SOI and B01 were based on the fit 
of the ISO SWS sp ectra of HM Sge o b tained in 10/01/1996 
and 05/16/1997 (|de Graauw et all 119961 : ISalama et all 
Il997fl. Assuming a period of the Mira variabl e of 527 days 
dWhitelockl . [1987) and using the ephemeris of lYudin et al.1 
(|l994h . they find that the observation were taken at phase 
0.65 (i.e. shortly after the photometric maximum) and 
at phase 1.08, close to the MIRA minimum, respectively. 
Both studies consider single and double shell models gen- 
erated with the DUSTY code. This public domain simula- 
tion code solves the problem of radiation transport in a cir- 
cumstellar dusty environment by integrating the radiative 
transfer equation in plane-parallel or s pheri cal geometries 
(jlvezic fc Elitzurl . 11996 al : llvezic et all . Il999f ). The studies 
consider different a priori (fixed) input parameters: 



5000L© and L/j Ot =9200L m luminosities con 



■ 



- S01 : L c 

sideri n g the period-luminosit y relation for Mira s by 
iFeastl i)l996l ) and the data of iMurset et all (|l997t ). re- 
spectively. They also take a syn t hetic stellar spectrum 
for the Mira from lLeieune et aT] (|l997f ) with a temper- 
ature of 3000K and loga=-0.29. 

— B01 : L COO ; = 10600L Q and L/j Ot =9000L Q luminosities 
considering the period-luminosity relation obtained 
from observations of Mira oxygen s tars in the Large 
M agellanic Cloud dFeast et al. . 119891 ) and the relation 
of lMurset k. Nussbaumerl (|1994 ). respectively. They im- 
pose that the shell outer boundary associated with the 
Mira and the WD is situated at thousand times the in- 
ner boundary radius. Finally, they take into account a 
2600K blackbody spectrum for the Mira. 

Both authors consider a 2xl0 5 K blackbody for the 
White Dwarf. In order to homogeneously compare the SEDs 
and visibilities of the models, we do not consider the syn- 
thetic stellar spectrum used by S01 and used a 3000K black- 
body stellar spectrum for the Mira. The standard size dis- 
tribution (MRN gr ain size distribution) as described by 
iMathis et all (|1977l) . n(a) oca with minimum and maxi- 
mum grain sizes of 0.005 /im and 0.25 /im respectively, and 
a chemical composition of warm silicate dust (W -Sil) with 
optical properties given bv lOssenkopf et all (|1992f ) . are also 
assumed by the 2 authors. 



3.2. Single shell models 

The SED of single shell models are shown in Fig. [3] The 
corresponding dust shell output parameters given by the 
DUSTY code for S01 and B01 are listed in Tab. 03 

Although both models provide good quality SED fits, 
their parameters are very different. One important diver- 
gence comes from the inner boundary temperature of the 
shell. S01 find a temperature (1600K) close to the conden- 
sation temperature of the grains while B01 find a weaker 
temperature (700K). The inner radius from B01 is 8 times 
further than the one found by SOI. The S01 visual optical 
depth is about twice larger (tv=29) than the one found by 
B01 (tv = 12.5) with a dust density power law coefficient of 
2.0 for B01 and 1.8 for S01. As a consequence, the S01 dust 
shell is more opaque than the B01 one, considering an iden- 
tical dust chemical composition (warm silicate dust). B01 
are satisfied with their single shell model whereas S01 think 
that the relatively high optical depth involved with their 
model ( >30 mag) is not compati ble with the JHK photom- 
etry by iKamath fc Ashokl (|1999( ) suggesting that the light 
from the Mira is not highly absorbed in near-IR. For this 
reason, the single shell model is rejected by S01, with the 
additional argument that such a model is probably simplis- 
tic in view of the disturbing presence of the hot companion. 



3.3. Double shell models 

Two shell models are an attempt to take into account the 
fact that each stellar component has its own effect on the 
dust generated in the Mira wind. The 2-component model 
can physically be justified with the scenario of a strong UV 
light source placed in the expanding shell of a red giant 
destroying the dust grains and generating an ionised nebu- 
losity. Radiation pressure from the companion will secondly 
slow the red giant wind and generate a compression wave 
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5 10 15 20 25 7 8 910111213 

wavelength (^m) wavelength (/inn) 



Fig. 3. Left: best fitting SOI single shell model (thick solid line) and best fitting B01 single shell model (thick dashed 
line), error bars correspond to the MIDI flux, both thin solid line are related to the ISO/SWS spectra. Right: close-up 
view of the best-fitted dust feature with the same labels. 

Table 2. Summary of the single shell DUSTY fit of HMSge data. Values marked with an asterix were pre-determined 
and fixed. 





Schild ct al. (2001) 


Bogdanov & Taranova (2001) 


Effective temperature (K) 


3000* 


2600* 


Luminosity (Lq) 


5000* 


10600* 


Distance (kpc) 


1.5 


2.6 


Central star diameter (mas) 


1.6 


1.8 


Grain chemical composition 


100% W-Sil* 


100% W-Sil* 


Grain size distribution 


MRN* 


MRN* 


Density power law coefficient 


1.8 


2.0* 


Inner boundary temperature (K) 


1600 


700 


Geometrical thickness (r^ n ) 


1000 


1000* 


Envelope inner radius (mas) 


3.5 


15.1 


Visual optical depth 


29 


12.5 


10/nn optical depth 


2.47 


1.77 



in the transition zone in which the conditions for dust for- 
mation are very favorable, so that the thin dust rich layer 
can be identified with its structure. For sake of simplicity, 
both studies assume that none of the component heats the 
dust shell of the other and that the total emitted flux is 
a simple addition of the individual fluxes from each shell. 
The best parameters found by S01 and B01 are gathered in 
Tab. [3 

S01 combine the two dust shells with a partition factor 
p such that the observed flux: 

Fobs = s(F\YDShell + pFMiraShell) (1) 

where the scaling factor s converts the model flux at the 
distance of HM Sge into the observed flux. One of the goals 
of the two shell modeling by S01 is to take into account 
the fact that the Mira in HM Sge is optically observable, 
at least in the near-IR. Hence, they fix an a-priori visual 
optical depth range between 0.5 and 3.5, much lower than 
in the single shell hypothesis. The authors also impose that 
the dust shell thickness associated with the White Dwarf 
should be geometrically thin (i.e. Y oll t/Yi ra <5) and there- 
fore could be approximated with a flat density distribution 
(x=0). The parameters of this geometrical configuration 
(Tab. [3]) consist in an optically thin, geometrically thick 
shell associated with the Mira, and a geometrically thin, op- 
tically thick shell associated with the White D warf. The S01 
2-component model is found satisfactory by ISchild et al.l 



(200lJ) but does not fit well shorter wavelengths dominated 
by the Mira spectrum. This can be attributed to the vari- 
ability of the Mira which is not taken into account by the 
model. Despite some discrepancies, the authors favor the 
2-component models approach. Fig. [4] shows their two shell 
best fitted model. The top-left part shows the S01 spectral 
energy distribution (SED) with the flux contribution of the 
individual shell, while its bottom-left part shows a close-up 
of the silicate dust peak around 10/im. 

The two shells model of B01 does not use any partition 
factor: it assumes that the flux from the central source is 
equal to the sum of two blackbodies located within a sphere 
of radius ri n . They argue that, given the small distance be- 
tween the components (i.e. 40 mas) compared to the size 
of their common dust envelope, this condition is approxi- 
matively satisfied. Contrary to the conclusion of S01, the 
results of the BOl's study given in Tab. [3] lead the authors 
to conclude that the 2 shell model is unable to describe 
the observed energy distributions at mid-IR wavelengths 
due to a twice larger squared deviation comparing to the 
single shell model. Moreover, the distance is a free param- 
eter of the model and B01 are lead to conclude that the 
overestimated distance values of the best fits (3.7kpc) is 
in contradiction with previous studies of HMSge. Fig. 2] 
shows their two shell best fitted model. The top-right part 
of Fig. [4] shows the B01 spectral energy distribution (SED) 
with the flux contribution of the individual shell, while its 
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Fig. 4. Top-left: best fitting SOI two shell model (thick solid line). Top-right: best fitting B01 two shell model (thick 
solid line), the thin dot-dashed line is the flux contribution of the Shell M (associated with the Mira), the thin dashed 
line is the flux contribution of the Shell WD (associated with the White Dwarf), error bars correspond to the MIDI flux 
and both thin solid line are related to the ISO/SWS spectra. Bottom: close-up view of the corresponding best-fitted dust 
feature with the same labels. 

Table 3. Summary of the two shell DUSTY fit of HMSge data. Shell "M" is the dust shell associated with the Mira 
component, Shell "WD" is the dust shell associated with the White Dwarf component. Values marked with a star were 
pre-determined and fixed. 





Schild ct al. (2001) 


Bogdanov & 


Taranova (2001) 




M 


WD 


M 


WD 


Effective temperature (K) 


3000* 


2xlO a * 


2600* 


2xl0 a * 


Luminosity (Lq) 


5000* 


9200* 


10600* 


9000* 


Distance (kpc) 


1.5 


1.5 


3.7 


3.7 


Central star diameter (mas) 


1.6 


5xl0 -4 


1.3 


2xl0" 4 


Grain chemical composition 


100% W-Sil* 


100% W-Sil* 


100% W-Sil* 


100% W-Sil* 


Grain size distribution 


MRN* 


MRN* 


MRN* 


MRN* 


Density power law coefficient 


1.6 


0* 


2* 


2* 


Inner boundary temperature (K) 


1400* 


800 


900 


900 


Geometrical thickness (n n ) 


1000 


1.1 


1000* 


1000* 


Envelope inner radius (mas) 


3.3 


128.7 


6.1 


40.0 


Visual optical depth 


1.5 


13 


8.8 


8.8 


10/im optical depth 


0.13 


1.11 


0.75 


0.75 


Relative contribution 


10 


1 


1 


1 



bottom-right part shows a close-up of the dust peak around 
10^m. 

We reproduce the mode ls of ISchild et al.l (|200ll ) and 
iBogdanov &: Taranoval (|200l[ ) . It is very interesting to point 
out that the values of the parameters and their respective 
conclusions are in complete opposition, B01 favor a single 
shell model while S01 favor a two shell model. The exten- 
sions of the dusty shells are very different from one model 



to another, and our mid-infrared interferometry data bring 
an important constraint that can help us to disentangle 
between them. 
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4. Visibilities of the models 

4.1. Description 

In this section, we use DUSTY to generate the synthetic 
spectrally-dispersed (7.5-13/zm) visibility profiles for the 2 
single shell and 2 double shell models of SOI and B01 for 
each projected baseline of the MIDI observations. Then we 
compare the synthetic visibilities with the mid-infrared vis- 
ibility measurements recorded. More details on the compu- 
tation of the synthetic visibilities with the DUSTY code 
can be found in Appendix A. 

4.2. Single shell visibility models 

Figure E] shows the SOI (solid line) and B01 (dashed line) 
single shell visibility models superimposed on the MIDI vis- 
ibilities (error bars) for each projected baseline. These plots 
show that, whatever the model, the visibilities are under- 
evaluated by comparison with the MIDI data. Nevertheless, 
the SOI model is relatively satisfactory, close to the ob- 
served visibilities, with a similar spectral shape. The B01 
model provides much lower visibilities. We point out that 
the visibility curve shape of the SOI model is purely related 
to the chemistry of the shell and not a consequence of the 
geometry of the source. The drop of visibility at lO^tm is the 
consequence of the silicate emissio n, leading to a more ex - 
tended object at these wavelengths ()Chesneau et all 12005). 
By contrast, the shape of the visibility curves from the B01 
model depends on the projected baseline with the wave- 
length of minimum visibility changing from 9.2 to 10/im. 
The shape of the B01 visibility curve reflects the change of 
angular diameter due to the silicates, but also the modu- 
lation of the visibility due to the ring resolved by the in- 
terferometer. The inner radius of the B01 model is large 
(angular diameter of 30 mas), and the dust shell is rather 
optically thin (rio Mm =1.8). Therefore, the dust emission is 
seen dominantly as a ring at low visibility level, whereas the 
SOI model is mostly seen as a compact source. Such mod- 
ulations from an ext ended dusty ring hav e already been 
observed with MIDI (jChesneau et all l2006t) . 

The shape of the observed visibility curves is relatively 
constant, because the shell is optically thick or because the 
inner rim radius is too small to be resolved by the inter- 
ferometer. The SOI model fulfills these two conditions: the 
10/im optical depth of the shell is 2.5, and the angular di- 
ameter of the inner rim is only 7 mas. 

Overall, the SOI model is impressively good, keeping 
in mind that the SED fits are usually not unique, and do 
not provide good predictions for the observed radii of dust 
emission. 

4.3. Double shell visibility models 

Computing the contribution from two separate shells and 
adding them to form a resulting SED with the DUSTY code 
is a simple operation provided that the fluxes are supposed 
to be independent. This strong hypothesis is not viable in 
the context of our interferometric observations that are also 
sensitive to the relative position of the two shells. Hence, 
it is necessary to add to the description of the two shells 
their projected separation and position angle of their cen- 
ters. Several options are possible: we position the centers 
of the two shells at the same position, or we position the 
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wavelength (/im) wavelength Qum) 

Fig. 5. Corresponding SOI (solid line) and B01 (dashed 
line) single shell visibility models superimposed on the 
MIDI visibilities (error bars) for the six projected baselines. 



centers at a given position and consider them as fixed pa- 
rameters, or we consider these news parameters as free and 
part of the fitting process. In order to limit the extent of 
the current study, we fixed the parame ters to the values of 
the binary found in lEvres et al.l (|2001h : the Mira and WD 
shells are positioned 40 mas apart at a PA of 130°. The 
detailed visibility computation is given in Appendix B. 

The direct consequence of the introduction of an ex- 
tended shell in the two components models is a significant 
decrease of the N band visibilities of the S01 and B01 model. 
In addition, some differences in the shape of the visibility 
curves can be noticed (see Fig|5]). The shapes of the visi- 
bilities of the B01 model exhibit a strong dependency on 
the position angle of the baseline. When the orientation 
of the projected baseline (~45°) is perpendicular to the 
binary separation, i.e. the projected angular separation be- 
tween the Mira and the WD is small (~3.5 mas for PA-46°) 
compared to the spatial resolution of the interferometer 
fFig lB.ip . the spectral shape and the level of the visibil- 
ities are naturally explained by the predominant correlated 
flux from the Mira shell seen by the interferometer. The 
unique effect from the extended (over-resolved at baselines 
larger than 30m) and optically thin (ri OMm =0.75) WD shell 
is to dilute the correlated flux to a low level with little spec- 
tral signature. When the orientation of the projected base- 
line gives a larger projected angular separation between the 
Mira and the WD (-23-35 mas at PA-75-100 ) compared 
to the spatial resolution of the interferometer (Fig lB.ip . 
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Fig. 6. Corresponding SOI (solid line) and B01 (dashed 
line) double shell visibility models superimposed on the 
MIDI visibilities (error bars) for the six projected baselines. 



the shift between the Mira and WD shells in the sky intro- 
duces an interferometric signal that mimics the one from 
a binary. Because the interferometer is able to resolve the 
2 shells with their own photocenter, a variation occurs on 
the slope of the visibility depending on the projected base- 
line. The binary effect is given by a sinusoidal dependency 
of the visibility shape versus the projected angular separa- 
tion of the 2 components. The visibility increases from 8 to 
13/im for the 15.8m baseline, then decreases for the 46.8m 
one, and increases again for the 59.3m one. This binary 
signal is not visible in the SOI model since the WD shell is 
over-resolved whatever the projected baseline and do not 
contribute to the correlated flux. The other differences be- 
tween the shape of the SOI and B01 visibility curves come 
from the flux balance between the emission of the two shells 
at various wavelengths, as can be seen in Figj4] The ratio 
of the two fluxes do not vary in the B01 model whereas 
a steep variation occurs between 9 and ll/mi in the SOI 
model. At 9.7/xm, the Mira shell dominates the flux leading 
to a local increase of the visibility. 

Clearly it is debatable whether these dust shell models 
associating two independent flux sources and dust distribu- 
tions are good approximations of complex sources as symbi- 
otic novae and provide a better way of understanding them. 
More complex geometries should be used, constrained by 
hydrodynamical models and high resolution imaging, as a 
further step for improving our knowledge of these sources. 



5. SOI single shell model improvement 

In this section we consider the SOI single shell model as the 
potentially best model able to account for spectra and vis- 
ibilities. We propose a pragmatic strategy of improvement 
that consists in small perturbations of the already satisfac- 
tory parameters of the SOI single shell model. We can see 
that the model is slightly too extended so that we need to 
spatially concentrate the emission of the dust shell. 

A solution consists in playing with some intrinsic pa- 
rameters of the dusty environment of HM Sge, i.e. the 
density distribution power law and the optical depth. 
Steepening the density power law (at constant mass) in- 
creases the dust density at the inner boundary with the 
density dropping faster. The emission comes from a shell lo- 
cated closer to the inner boundary and, consequently, less 
resolved by the interferometer. A decrease of the optical 
depth (by decreasing the total amount of mass) leads to a 
more transparent dusty environment, so that deep seated 
layers contribute significantly to the correlated flux seen 
by the interferometer. These two changes do not affect the 
SED in the same way. The dust density distribution can 
only be constrained from the SED beyond 10^m since the 
matter distribution mainly determines the slope of the re - 
distributed energy by dust emission (|Lobel et all I1999D . 
On the other hand, a variation of the optical depth affects 
the whole shape of the curve, in particular the silicate fea- 
ture. However, playing with both parameters, we find that 
the dust density distribution law has much more influence 
on the visibilities than a significant decrease of the optical 
depth. Moreover a decrease of the optical depth changes 
the whole SED curve that quickly does not fit the ISO and 
MIDI fluxes anymore. Diminishing the density power law 
coefficient increases the visibility but differences between 
the flux data and the model arise at longer wavelengths, 
with a too small flux of the model compared to the ISO 
data. It could be argued that this is due to the appearance 
of a beam size effect toward wavelengths larger than 15/im, 
the ISO telescope gathering light from a large field-of-view. 
But the MIDI visibilities suggest a compact dusty environ- 
ment that easily fits into the 8m telescope beam of 300 mas 
and the ISO emission is most probably also spatially com- 
pact at larger wavelengths. 

The perturbed model has a density power law coefficient 
of 2.2 instead of 1.8 and a 10/j,m optical depth of 2.15 in- 
stead of 2.47 with the remaining parameters equal to those 
of the SOI single shell (see Tab. [5]). The fits are rather good 
but the optical depth of the perturbed model remains high. 
This model is subject to the same critics than for the SOI 
one: the shell prevents the flux from the Mira to be really 
detectable in near-IR. We tested a last change in order to 
improve the quality of the fits between 8 and 10/iin. One 
can notice that the model visibility slope between 8 and 
10//m is too strong revealing a too much extended struc- 
ture at these wavelengths. Introducing a contribution from 
an optically thinner material than silicate such as carbon 
may allow to decrease the size of the object at these wave- 
lengths. 

Carbon in the form of graphite has b een often involved 
to co nstruct the early SED models (jKwok &: Purtonl . 
I1979D . and its presence was justified by the fact that this 
material would come from the White Dwarf ejecta dur- 
ing the outburst of 1975 or due to particular dust nucle- 
ation processes that occur under the non-equilibrium phys- 
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Fig. 8. Top: Best perturbed visibility profile model super- 
imposed on the MIDI visibilities at 8.5, 10 and 13/im (from 
the left to the right). Bottom: Corresponding visibility dif- 
ference between the model and the MIDI data. 



ical conditions that prevail around the WD. An other ar- 
gument is also provided bv lBrvan fc Kwo"kl (|1991l ) in their 
analysis of the time variable 0.22/im feature in IUE spec- 
tra that they attribu t e to c arbon-based material (see also 
iMunari fc Whitelockl (|1989l )). If present, this carbon con- 
tamination must be small and we developed a test model in 
which the dust composition is a mixture of 90% of warm sil- 
icate s (lOssenkopf et all I1992D and 10% of amorphous car- 
bon (|Hannerl . ll988l ). Fig. shows the best perturbed S01 
single shell model with the new mixture: the density power 
law coefficient is 1.8 and the 10/zm optical depth has de- 
creased to 1.89. The other parameters of this model keep 
the same values of the S01 model contained in Tab.[H Fig. [5] 
shows the best perturbed visibility profile model superim- 
posed on the MIDI visibilities at 8.5, 10 and 13/tm and the 
visibility difference between the model and the MIDI data. 
Fig. [S] shows the corresponding best perturbed intensity 
profile model at the same wavelengths. 



6. Discussion 

6.1. Spherical models: single versus double shells 

The main result of the MIDI observations is that the dust 
is spatially concentrated around the Mira star and that the 
MIDI single-dish fluxes is directly comparable to the ISO 
one. Using a Gaussian model, we estimated the HWHM of 
the dust in the 8-9/im and 11-12/im range to 7.8±1.3mas 
(12AU) and 11.9±1.3mas (18AU) respectively. 

We compared the MIDI data with the predictions of var- 
ious single and double spherical shells spectrophotometric 
models and favored the single shell models as they pro- 
vided the most spatially concentrated emission. Gaussian 
model fitting of our best model provide the following Half 
Width at Half Maximum: 7.2 mas at 8.5/im, 13.2 mas at 
10/im, and 11 mas at 13/im (see Fig. [9]). These extensions 
are much smaller than the separation between the Mira and 
the W D of 40 mas (60 AU ) determined by the HST observa- 
tions (jEvres et all l200"Th . At 20 mas, the local intensity of 
the envelope is less than 3% than the one at 4 mas, and the 
amount of dust that is located in the vicinity of the WD is 
limited. 

The mass- loss rates of the best S01 and the per- 
turbed models are 7.7xl0~ 6 and 7x 10~ 6 M w /yr, respec- 
tively. Similar estimations can be found in [Richards et al.l 




Fig. 9. Top: Best perturbed intensity profile model at 8.5 
(solid line), 10 (dotted line), and 13/im (dashed line). The 
vertical lines indicate the HWHM (at 8.5 (solid line), 10 
(dotted line) and 13/im (dashed line)) of the best per- 
turbed intensity profile using a Gaussian fitting. Bottom: 
The mid-infrared images (at 8.5 (left), 10 (middle), and 
13/im (right)) predicted by the best perturbed model given 
in a logarithmic scale. 



(1999]). lNussbaumer fc Vogell (|1990t ). and lKenvon fc Garcial 
(|1989fi . 

The difficulty encountered in the fitting process is that 
a large amount of dust is required to explain the 10-100/im 
flux and that the Mira photosphere should be visible in the 
near-IR. The solution proposed by (|Brvan fc Kwokl Il991h 
is to use a large inner radius, 10. 5R* in the case of HM Sge, 
but such a large inner radius would increase strongly the 
mid-IR extension of the system and is not compatible with 
the data. In the other hand, the best perturbed S01 model 
is able to fit many observations in first order but is still far 
from providing a satisfactory fit of the near-IR and long 
wavelength SED and is also too optically thick. 

At this point, the main assumption of spherical symme- 
try has probably reached its limits. Although there is con- 
siderable evidence for departures from spherical symmetry 
at large scales (>1"), the primary constraints on the mod- 
els from the MIDI mid-IR visibilities and spectra are not 
sufficient to warrant more complex models at small scales. 
T he use of more com plex models, like the one developed 
in lLopez et alj (|l997l ) would help us to better constrain the 
dust density distribution of the system. However, such mod- 
els have to be well defined through a large uv coverage and 
require a good knowledge of the wind-wind collision zone. 
We clearly lack this information with our limited data set 
and we consider this modeling out of the scope of the pa- 
per. Nevertheless, we investigate further whether any sign 
of departure from spherical symmetry of the mid-IR source 
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Fig. 7. Top-left: best fitting of the perturbed SOI single shell model (thick solid line), error bars correspond to the MIDI 
flux and both thin solid line are related to the ISO/SWS spectra. Top-right: close-up view of the best-fitted dust feature 
with the same labels. Bottom: Corresponding single shell visibility model superimposed on the MIDI visibilities (error 
bars) for the six projected baselines. 



is seen by comparing the circumstellar dust extension ac- 
cording to the baseline orientation. 



6.2. Departure from spherical symmetry 

Our projected baselines cover position angles from 42° to 
105°, allowing us to put some constraints on the degree 
of spherical geometry of the dust. It is also of prime im- 
portance to determine as best as possible whether a binary 
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signal is present in the MIDI data and whether the position 
angles of the projected baselines are aligned with a partic- 
ular direction reported for the HM Sge binary system from 
radio, optical or polarization data. 

Let us recall first that even tho ugh the WD component 
has a high luminosity (~10 000L^) (jMurset fc Nussbaumerl . 
1 1994 iMurset et all Il997t) and the binary system separa- 
tion (~40mas) is optimum for a detection with the MIDI 
instrument, its low mid-IR flux prevents MIDI to directly 
detect its signal (F^/F^^IO 5 ). The mid-infrared 
flux coming from an accretion disk would add to the WD 
one, but the high temperature of such a structure makes 
it also too faint in mid-IR. The WD and a putative accre- 
tion disk could potentially perturb the dusty environment 
leading to a large scale signature mimicking a binary sig- 
nal. Such a signal is probably observed in the case of R Aqr 
(jTuthill et all I2000T ). that the authors attribute to a com- 
pact mid-IR source representing 3% of the N band flujjf] 
at 700 mas from the cool star (this separation would scale 
down to below 100 mas at 1.5kpc). Such a feature might 
be some localized concentration of dust, a local warming or 
a global asymmetry of the dust shell as also seen and dis- 
cussed in Secf4]in the frame of the B01 double shell model 
that predicts such a binary signature. 

As previously said, the MIDI visibilities are rather 
s mooth and no hin t of su ch oscillation (as seen for instance 
in IChesneau et all (|2006l )) is found in the data (the shape 
of the visibility profile may uniquely be related to the dust 
shell). It is possible to be more sensitive to a weak signal of 
this kind by using the differential phase. Fig. fTOl shows the 
calibrated differential phase of the object evaluated from 
all the determined projected baselines (dotted lines). Error 
bars correspond to the rms of the differential phase over 
all the projected baselines. The rms values of the differen- 
tial phase (±4.26° at 8.5/im, ±1.79° at 10.5/mi and ±1.66° 
at 12.5/Ltm) indicates that there is no obvious signature of 
any astrophysical effects. The most discrepant curves are 
not correlated to a particular baseline length nor position 
angle and seem to reflect fast atmospheric changes during 
observations that were not fully corrected by the calibrator. 

In order to estimate further the constraints provided 
by such measurements, we performed a perturbation of our 
best single shell model by a point-like structure at vary- 
ing flux contrast, distance and angle. Our tests allow us to 
discard any compact mid-IR source representing more than 
about 1% of the total flux (i.e. 0.5 to 1 Jy between 8/j.m and 
13/im) at a distance ranging from 20mas to 80mas from the 
Mira staiQ. 

Let us now consider the level of large scale asymme- 
try for the source. Fig. QT] shows the equivalent HWHM 
of the Gaussian intensity distribution according to the ori- 
entation of the projected baseline on the sky (UT2-UT3 
and UT3-UT4). The figure is centered on the Mira and 
present s the binary a xis di re ction (around 1 20-13 0° accord- 
ing to lEvres et alj (|200lh : ISchmid et all (120001)1 located 
beyon d the window at 40 mas north-west (jEvres et all 
2001). The MIDI baselines at PA~45° are oriented al- 
most perpendicular to binary separation, in t he direction 
of the outflows evidenced by various authors ((Sol . fl98l 



Note that this signal is also seen in recent MIDI data. 
4 With a separation of 40mas, at full sinusoidal pattern should 
be seen in the spectrally-dispersed MIDI visibilities and differ- 
ential phases from 8(j,m to 13/xm. 
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Fig. 10. Calibrated differential phase of the object evalu- 
ated from all the determined projected baselines (dotted 
lines). Error bars correspond to the rms of the differential 
phase averaged over all the projected baselines. 



Hac k fc Paresce L 19931: Richards et all Il999t ICorradi et all 
11999 ; Eyres et aLL 200l[ ). On the contrary, the baselines at 
PA^100° are favorably oriented to probe the largest exten- 
sion of the dust, if this dust is elongated in the preferential 
direction of the shadow of the Mira, that should efficiently 
protect the dust from the hard emission of the WD. The 
orbital period of HM Sge is unknown but certainly rather 
long and we assume that it is barely unchanged since 1998, 
date of the latest polarization and HST observation^. 

Using the Gaussian HWHM estimations, the following 
ratio between the HWHM at 105° and the ones at 42° can 
be computed at 8, 9, 10 and 13/xm, respectively: 0.72 



0.80 



+0.20 



+0.32 
-0.19' 



, 0.84t°;^ and 0.91±g;J|. The extension of the shell 
is larger in the perpendicular direction than in the binary 
axis. The level of discrepancy is larger at 8/zm and steadily 
decreases toward longer wavelengths, with HWHM from 5 
to 10 mas at 8/xm, to about 11 to 13 mas at 13/im. 

The wind from the cool star is likely to have an equa- 
torial density enhancement. Spectro-polarimetric observa- 
tions show that continuum from th e Mira star is intrinsi- 
cally polarized (|Schmid et all l2000f ) and this polarization 
would be, according to the author, most likely due to light 
scattering in an aspherical photosphere or wind. 
Several geometries can be envisaged: 

— A disk-like geometry seen almost edge-on. We expect 
to detect a strong difference between the MIDI visibil- 
ities with perpendicular baselines that is not observed. 
Moreover, depending on the inclination and the disk 
vertical extension, a large circumstellar absorption of 
the Mira light is expected. 

— A disk-like geometry seen almost pole-one offer an ele- 
gant solution to the opacity problem: the light from the 
Mira star can reach us through an optically thin path 
whereas a dense and compact disk seen pole-on emits a 



5 We notice that the proper mo t ions reported for some ex- 
tended features in Richards et alj (1 19991 ) is typically 2.5° per 
year from north to east. If this estimation can be extrapolated 
to the binary motion, then the rough expected position angle 
of the projected binary plane should be around 145-150° at the 
time of the MIDI observations. 
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large infrared excess. This geometry is relatively sym- 
metrical and exhibit similar extensions whatever the di- 
rection of the baseline is. 

A slightly more complex geometry in which the 
Mira is wrapped into a sh ock cone (Fig[T2"l) as pre - 



1995) 



2005). 



sented in several studies Fbrmiggini et al.J 
iMastrodemos fc Morris! (|l999f ): lKennv fc Tavlorl 
The density is not uniform and the envelop may be thin- 
ner (and patchy) in the direction perpendicular to the 
binary direction, allowing the Mira light to reach us, 
and a large dust content, mostly symmetrical to reside 
close to the Mira. This geometry is identical for pole-on 
or edge-on inclinations. 

The probably high degree of symmetry of the source 
favors the two last scenario in which most of the dust resides 
in a spherical/disk- like pole-on envelope around the Mira, 
that is currently only slightly perturbed by the current hard 
flux of the WD. Such a low level of perturbation may appear 
surprising considering the strong uv source around the cool 
star, but the separation is rather large (~60AUs), and the 
shielding of the uv flux by the dense and slow Mira wind is 
probably im portant (see the case of the symbiotic nebula 
PN M2-9 in lLivio fc Soken pOOlh ). 

A continuum of inclinations is possible but the observa- 
tions of a symmetrical optical and radio outflow suggest a 
high inclination. In the frame of the wind-wind collision ge- 
ometry, the larger extensions perpendicular to the binary 
axes could be explained by a somewhat opened aperture 
angle of the shock cone. At 8/xm, hot dust emits preferen- 
tially close to the inner diameter of the Mira shell, but an 
additional emission component from the shock front may 
provide the asymmetry observed. As the contribution of 
the shock component decreases from 8 to 13/xm, the bulk 
of the emission is more centered around the Mira and the 
level of asymmetry decreases (see FigfT^]). 

We note that the wind-wind collision hypothesis has 
been investigated and model e d many times in the litera- 
ture of HMSge (lEvres et all. l200ll; [Richards et all Il999t 
iGirard fc Willsonl . 119871 : IWillson et al.l . ll984allbD . The tech- 
nique of mid-IR interferometry is a good opportunity to 
extensively study such a structure, but the present paper 
shows that two-telescopes observations, as currently per- 
formed by MIDI, are not sufficient to constrain its shape 
efficiently. What is clearly needed is an interferometric im- 
ager able to recombine simultaneously the light from the 
4 UT or AT telescopes in the mid-IR. Such a recombiner 
called MATISSE is propose d for the second ge neration in- 
strumentation of the VLTI (|Lopez et all l2006h . 



7. Summary and conclusion 

We have presented new high spatial resolution observa- 
tions of the dusty environment of HMSge, consisting on 
MIDI/ VLTI spectra and visibilities. Using the strong con- 
straints provided by optical interferometry, we reassessed 
previous interpretations based only on the modeling of 
the spectral density distribution. In the case of the com- 
plex symbiotic systems, the spe ctral energy densit y can 
be modeled in different wa ys . ISchild et all (1200 ll ) and 
Bogdano v fc Taranoval (|2001[ ) have worked the same year 
on the ISO/SWS spectra of HM Sge. Both take into account 
single and double shell models, as an attempt to take into 
account the spatial complexity of the source. Both used 




Fig. 11. Equivalent HWHM of the Gaussian intensity dis- 
tribution (in mas) dependence on the projected baseline 
position angle on the sky. These HWHM are taken at 8 
(plus), 9 (cross), 10 (triangle) et 13/um (square). The figure 
is arbitrarily centered on the Mira and presents the binary 
axis (dashed line) (PA=130°) (lEvres et al.l.l2001l) where th e 
position of the White Dwarf at 40 mas (|Evres et al.1 . l200l[ ). 
represented by the arrow, is out of the window. The thick 
line, perpendicularly to the binary axis, rep resents the po- 
sition angle of the Raman line polarization ( Sc hmid et all 
2000). The black filled circle represents the dimension of 
the Mira determined with the best perturbed model. The 
dotted circle represents the dust shell inner boundary di- 
mension determined with the best perturbed model. The 
two 8 and 13^tm solid ellipses have a minor axis of 5.6 and 
11.1 mas, and a major axis of 7.8 and 12.5 mas, respectively. 



the same radiative transfer code DUSTY. However, they 
do not converge to the same value of parameters and their 
respective conclusion are in complete opposition. The inter- 
ferometric study of this object shows that the morphology 
of the dust was cl ose enough to the single shell model of 
ISchild et al.l {2001) considering that the circumstellar envi- 
ronment of the symbiotic system was compact and optically 
thick (an angular inner radius of the dust shell located at 
3.5 mas and a visual optical depth of 29). 

An improve ment of the model c onsisting in small per- 
turbations of the lSchild et ail (|200lD parameters shows that 
a good fit of both SEDs and visibility curves is obtained 
steepening the density power law coefficient from 1.8 to 
2.2 and decreasing the 10/xm optical depth from 2.47 to 
2.15. The change of these two parameters gives a larger 
transparency of the circumstellar medium and allows the 
emission to come from deeper layers with the aim to spa- 
tially concentrate the emission of the dust shell. Finally 
the insertion of carbon material in the dusty layers allows 
to diminish a too much extended structure between 8 and 
10/im. Such a chemical component could come from the 
White Dwarf ejecta during the outburst of 1975 or from 
particular dust nucleation processes that occur under the 
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Fig. 12. Schematic representation of the symbiotic system 
in which the Mira (in the center of the figure) is wrapped 
into a shock cone. In the top on the right is the White 
Dwarf. The arrows surrounding the 2 components indicate 
the 2 winds and the parallel arrows symbolize the radiation 
reaching the shock front component (gray filled zone). The 
contribution of the shock component decreases from 8 to 
13/im simultaneously with the level of asymmetry. The 3 
circular regions correspond to the dust shell extension at 8 
(solid line), 10 (dotted line), and 13/xm (dashed line). 



non-equilibrium physical conditions that prevail around the 
White Dwarf. 

The compact model was rejected bv lSchild et al.l (|200ll ) 
because of its high optical depth not compatible with non- 
absorbed emission of the Mira in the near-IR. Our best 
spherical model has to face the same objection, but this is- 
sue could be solved assuming a more complex spatial distri- 
bution of the dust. Our MIDI data allowed us to investigate 
the departure of spherical symmetry using different base- 
line orientations. The extension is stronger in the direction 
perpendicular to the binary axis. This asymmetry is largest 
at 8/im and steadily decreases toward longer wavelengths. 
We speculated on the origin of the departure. We propose 
that the Mira is wrapped into a shock cone where the den- 
sity of the envelope is thinner and patchy in the direction 
perpendicular to the binary axis. This could allow the Mira 
light to reach us, and a large dust content, mostly sym- 
metrical to reside close to the Mira. Further VLTI/MIDI 
observations at 15 different projected baselines with ATs 
are planned to confirm this geometrical scenario. 

The MIDI observations must be included in the frame 
of the long-term evolution of the dusty envelope of HM Sge 
associated with the variability of the conditions of dust sur- 
vival in the system, and hence the luminosity and temper- 
ature of the White Dwarf, and also the relative position 
of the Mira and the WD in the orbit. The current dusty 
envelope from the Mira appears only slightly perturbed by 
the hot companion. As the time passes, the impact of the 
WD ionizing flux on the Mira dust shell should decrease 
further for this large separation binary system unless that 



in the course of the putatively eccentric orbit a decrease of 
the separation occurs, or that a large accretion disk builds 
up. 
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Fig. A.l. Geometry used for the determination of the 
emergent intensity from the shell outer boundary (P) and 
received under the angle e by the observer located at O, 
distant of D from the center C of the central source. The 
point r is the impact parameter; ri n is the inner bound- 
ary radius, and £i n is the angular radius of the inner shell 
boundary. 



where £® n is the angular radius of the inner boundary when 
Fbol = 10~ 6 W.m~ 2 and Ftoi is the observed bolometric 
flux of the source. 

Therefore, the spatial frequency / is given by 



l F p 

r bol 



2e* 



2f l 



Ft 



hoi 



(A.6) 



DUSTY can produce maps of the visibility at up to 20 
wavelengths. The visibilities of 20 spectral channels arc gen- 
erated between 7.5 and 13/im corresponding approximately 
to the number of independent channel provided by MIDI 
with the low resolution mode (prism, R^20). 



Appendix A: Determination of the model visibilities Appendix B: The binary separation 



The DUSTY code is used to get a theoretical model of the 
normalized radial intensity distribution at specified wave- 
lengths. For each model, DUSTY provides the radial in- 
tensity (in Jy.arcsec -2 ) versus the invariant dimcnsionless 
reduced impact parameter p defined by 



r 

P = — 

Tin 



(A.l) 



In the frame of a two shell model, the relative position of 
the centers has to be taken into account in the computation 
of the visibilities. We adopted a binary separation and a 
position angle of the binary axis es timated of 40 mas and 
130° respectively (|Evres et all 12001). Putting the Mira star 
at the center of reference (see Fig lB.ll) , the expression of the 
visibility can be written as: 



where r is the impact parameter and r„ is the radius of 
the inner boundary (see Fig lA.ip . 

Taking into account that the star in the code is nor- 
malized at luminosity = 10 4 Lq, the radius Ti n is given 

by 



(A.2) 



where rf n is the radius of the inner boundary when = 
1O 4 L and L* is the luminosity of the star. 

Therefore, the impact parameter r is given by 





(A.3) 



where the luminosity corresponds to the one given by the 
two authors for the Mira and the White Dwarf (see Tab(2] 
and Tab ED. 

DUSTY also allows us to generate the theoretical vis- 
ibility curves by calculating the Hankel transform of the 
circularly symmetric intensity distribution. The code gives 
the visibility versus the invariant dimensionless scaled spa- 
tial frequency 4> defined by 



= /.2e„ 



(A.4) 



where / is the spatial frequency (in arcsec 1 ) and Si„ is the 
angular radius of the shell inner boundary (see Fig |A.l[) . 



Taking into account that DUSTY considers a star with 
an observed bolometric flux = 10~ 6 W.m~ 2 , the angu- 
lar radius of the inner boundary is given by 



' Fboi 

r bol 



(A.5) 



V = 



totCi 



(B.l) 



y/ch + C wd + 2CmCwdcos(2tt^ Po cos{0i - )) 

Fm + F\YD 

where Ctot is the total correlated flux of the Mira plus the 
White Dwarf (C tot = C M + C WD e~ 2m ^- poCOS ^'- eo '>) and 
C£ ot is its conjugate. B=(Br,6i) are the coordinates asso- 
ciated to the interferometer projected baseline on the sky 
where Bj is the baselength and 0i its position angle from 
the north to the east. ~f?=(po;8o) are the coordinates as- 
sociated to the binary components where po is the binary 
separation and 9q the position angle of the binary axis from 
the north to the east. Finally, F to t represents the total flux 
of the Mira plus the White Dwarf (Ftot = F M + F WD ). 
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Fig. B.l. Schematic view of the separation between the bi- 
nary components. po=40mas is the distance between the 
Mira and the White Dwarf (W D) and fl =130° is the po- 
sition angle of the binary axis (|Evres et all . |2001| ). 7 are 
the position angles (at 45° and 100° from the north to the 
east) of the interferometric baselines from the sky. The thick 
lines correspond to the projected separations of the binary 
components on the sky. 
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